Draft version October 17, 2012 

Preprint typeset using I^T^^X style emulatcapj v. 5/2/11 



COLLISIONLESS SHOCKS IN A PARTIALLY IONIZED MEDIUM. 
II. BALMER EMISSION 

G. MoRLiNo\ R. Bandiera\ p. Blasi^^, E. Amato^ 

^INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5, 50125 Firenze, Italy and 
^INFN, Laboratori Nazionali del Gran Sasso, Assergi (AQ), Italy 
Draft version October 17, 2012 

ABSTRACT 

Strong shocks propagating into a partially ionized medium arc often associated with optical Balmer 
lines. This emission is due to impact excitation of neutral hydrogen by hot protons and electrons in 
the sh ocked gas. The st ructure of such Balmer-dominated shocks has been computed in a previous 
paper (jBlasi et al.l 12012}) . where the distribution function of neutral particles was derived from the 
appropriate Boltzmann equation including coupling with ions and electrons through charge exchange 
and ionization. This calculation showed how the presence of neutrals can significantly modify the shock 
structure through the formation of a neutral-induced precursor ahead of the shock. Here we follow up 
on our previous work and investigate the properties of the resulting Balmer emission, with the aim of 
using the observed radiation as a diagnostic tool for shock parameters. Our main focus is on Supernova 
Remnant shocks, and we find that, for typical parameters, the Ha emission typically has a three- 
component spectral profile, where: 1) a narrow component originates from upstream cold hydrogen 
atoms, 2) a broad component comes from hydrogen atoms that have undergone charge exchange 
with shocked protons downstream of the shock, and 3) an intermediate component is due to hydrogen 
atoms that have undergone charge exchange with warm protons in the neutral-induced precursor. The 
relative importance of these three components depends on the shock velocity, on the original degree 
of ionization and on the electron-ion temperature equilibration level. The intermediate component, 
which is the main signature of the presence of a neutral-induced precursor, becomes negligible for 
shock velocities ^ 1500 km/s. The width of the intermediate line reflects the temperature in the 
precursor, while the width of the narrow one is left unaltered by the precursor. In addition, we show 
that the profiles of both the intermediate and broad components generally depart from a thermal 
distribution, as a consequence of the non equilibrium distribution of neutral hydrogen. Finally, we 
show that a significant amount of Balmer emission can be produced in the precursor region if efficient 
electron heating takes place. 

Subject headings: acceleration of particles - atomic processes - lineiprofiles - ISM: supernova remnants 



1. INTRODUCTION 

It has been shown (jChevalier fc Ravmondl [T978D that 
optical spectra dominated by Ha and other Balmer 
lines, as observed in some historical Supernova Remnants 
(SNRs), may arise when an astrophysical shock propa- 
gates through a partially ionized medium. 

The Balmer emission, observed from these so-called 
Balmed-dominated shocks, provides a powerful diagnos- 
tic tool to investigate the conditions existing in the shock 
vicinity. The Ha lines typically show two components, 
resulting from excitation of neutral hydrogen due to the 
interaction with hot protons and electrons in the shocked 
gas: a narrow-line component, whose width is character- 
istic of the cold interstellar medium and that has been 
explained as the result of direct excitation of neutral hy- 
drogen atoms; and a much broader component, associ- 
ated to a second population of hydrogen atoms, created 
by charge-exchange (CE) processes between the cold, still 
unshocked hydrogen and the shocked protons. These hot 
atoms can be produced in an excited state or can be ex- 
cited by subsequent collisions with protons or electrons. 
Hence the line width of the broad component traces the 
thermal velocity of shocked protons and can be used to 
infer the shock velocity. Combining this estimate with 
proper motion measurements, one can estimate the dis- 



tance to the object. 

Besides the shock speed, Balmer lines also represent 
a unique tool to investigate the plasma physics of coUi- 
sionless shocks. If both narrow and broad components 
are detected, the relative intensity of the two lines can 
be used to infer the ratio of electron-to-ion tempera- 
ture just behind the shock, providing i nformation on the 
electron-ion equilibration mecha nisms (jChavamian et aP 
[2007HvMrAdelsberg et an[200l . 

One of the most intriguing aspects of Balmer emission 
is related to the possibility of using the line shape and 
its spatial profile to check the efficiency of SNR shocks in 
accelerating cosmic rays (CRs). If CR acceleration is tak- 
ing place in an efficient way, then the widths of both the 
narrow and broad line may be affected. In fact, when a 
sizable fraction of the ram pressure is channeled into non- 
thermal particles, the plasma temperature behind the 
shock is expected to be lower, and this should reflect in a 
narrower width of the broad Ha line. On the other hand 
efficient particle acceleration also leads to the formation 
of a CR-induced precursor upstream, which heats the 
ionized plasma before the shock. If the precursor is large 
enough, CE can occur upstream leading to a broader nar- 
row Balmer line. Remarkably, both signatures seem to 
have been observed i n Bal mer-dominated shocks. For ex- 
ample [HeldeZer^ (|2009f ) combined proper motion mea- 
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surements of the shock and broad Ha hne width for the 
remnant RCW 86 to demonstrate that the temperature 
behind the shock is too fow, thereby concluding that a 
sizable fraction of the energy is being channelled into cos- 
mic rays. However, qualitatively similar features can also 
arise from different physical processes. Therefore their 
observation can only be turned into reliable information 
on the shock properties after a quantitative physical de- 
scription of the phenomenon is provided. 

The basic theory of collisionless shocks in the 
pre sence of neutral par t icles was first developed 
bv iChevalier fc Ravmondl (119 781) and iChevalier et aD 
(I1980D and furth er refi'^d bv iSmith et al.l (|199lD "a^ 
iGhavamian et all ()200lD . These papers were however 
characterized by similar rather important limitations, 
mainly due to the assumption that the distribution func- 
tions of both populations of neutrals are Maxwellians. 
In fact, the main difficulty in describing the structure 
of a collisionless shock propagating in a partially ion- 
ized medium is that neutrals have no time to reach ther- 
malization and cannot be treated as a fluid. Steps for- 
ward in rel axing the fluid assumption hav e been made by 
iHeng fc Mc Crav (2007) and van Adelsb erg et al.l (|2008f ) 
which considered the effect of multiple CE events on the 
distribution function of hot neutrals. However these au- 
thors limit their calculations to the region downstream 
of the shock, and consider only the volume integrated 
distributions. 

A proper description of the effects of the interactions 
between neutrals and ions on the neutral distribution 
function requires a fully ki netic approach. Th is approach 
has been implemented by iBlasi et "all ()2012|) (hereafter. 
Paper I) where we derived simultaneously the neutral 
distribution function and the hydrodynamic quantities 
for ions (which are, instead, treated as a fluid), both up- 
stream and downstream of the shock, by solving a Boltz- 
mann equation including CE and ionization terms. The 
main result of Paper I is that of providing a mathemat- 
ical and physical description of what we call the neutral 
return flux: when fast, cold neutrals undergo CE interac- 
tions with the slower hot ions downstream of the shock, 
some fraction of the resulting hot neutrals can cross the 
shock and move upstream. The relative velocity between 
these hot neutrals and the upstream ions triggers the on- 
set of CE and ionization interactions that lead to the 
heating and slowing down of the ionized component of 
the upstream fluid. The system then tends to develop a 
shock precursor, in which the fluid velocity gradually de- 
creases from its value at upstream infinity, and even more 
important, the temperature of ions in the upstream re- 
gion increases as a result of the energy and momentum 
deposition of returning neutrals. 

The existence of a neutral return flux from downstream 
to up strea m was previously in entioned by ISmith et al.l 
(fT99l and iHester et all (|l99l as a possible way to ex- 
plain the anomalous width of narrow Balmer line s ob- 
served in some SNRs (see e.g. lSollerman et al]|2003[ ). but 
no explicit calculation was carried out. 

The width of such lines is in the 30-50 km/s range, im- 
plying a pre-shock temperature around 25000-50000 K. If 
this were the ISM equilibrium temperature there would 
be no atomic hydrogen, implying that the pre-shock hy- 
drogen is heated by some form of shock precursor in a 
region that is sufficiently thin so that coUisional ioniza- 



tion equilibrium cannot be established before the shock. 
Several mechanisms have been pro posed to ex plain the 
broadening of the narrow line (see Ee^[200l for a re- 
view), but most of them can be ruled out on theoretical 
grounds, leaving a CR precursor and/or a neutral pre- 
cursor as the most probable origin. 

A first attempt at investigating the broadening of the 
narrow line comp onent induced by th e neutral precur- 
sor was made by iLim fc Ragal ()1996[) . using a simpli- 
fied Boltzmann equation for neutrals in one dimension in 
both physical and velocity space. The calculations were 
carried out for a shock speed of 225 km/s, and varying 
the initial ionization fraction from 0.5 to 0.99. The nar- 
row line width does not show any appreciable change due 
to the return flux. 

In the present work we use the theory developed in 
Paper I and calculate the profile of Balmer line emission 
in the presence of neutral return flux. In particular we 
show that the neutral return flux is responsible for the 
emission of a third intermediate line component, in ad- 
dition to the narrow and broad ones. This intermediate 
line is produced by hydrogen atoms that have undergone 
CE with warm protons in the neutral precursor. Inter- 
estingly, there are observations which suggest the exis- 
tence of such intermediate component (jChavamian et ahl 
l2000f ). although the results might also be due to projec- 
tion effects since the emission region is morphologically 
rather complex. At the present time, spectral and spatial 
information are not sufficient to disentangle the physical 
effect that we describe here from geometrical and pro- 
jection effects which could modify the line profiles. We 
also show that the neutral-induced precursor is not able 
to broaden th e width of narrow H a line, confirming the 
first finding of ILim fc Ragal ((19961) . 

The paper is organized as follows: in 121 we summarize 
the kinetic approach developed in Paper I and use it to 
describe the shock structure in the presence of neutrals. 
We improve on previous work by both adding electrons in 
the shock dynamics and including their contribution to 
the ionization of neutrals. In SjSjwc write down the basic 
equations for the calculation of the Balmer line emission 
and in i|4| wc illustrate the main results in terms of spa- 
tial profiles of the total emission, of line profile, as well as 
line intensity ratios. We also compare theoretical profiles 
with profiles obtained by simulated observations, in or- 
der to derive the observational requirements necessary in 
order to be able to detect the intermediate component, 
and deviations from gaussianity in general. We conclude 
in gH 

2. PHYSICAL MODEL 

The basic model wc consider consists of a plane-parallel 
shock with velocity Vsh that propagates into a partially 
ionized proton-electron plasma along the z direction, 
with a given fraction of neutral hydrogen at upstream 
infinity. We neglect the presence of helium and other 
heavier chemical elements. 

Protons and electrons are assumed to behave as fluids 
with temperatures Ti(z) and Te{z), respectively, with the 
same bulk velocity, Vi(z) = Ve{z), and the same num- 
ber density, ni(z) = ne{z). Their distribution functions, 
fi{v,z) and fe(v,z), are assumed to be MaxwcUian at 
each position z. Neutral hydrogen interacts with protons 
and electrons through CE and ionization. The hydrogen 
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distribution function, /^{v^z)^ can be described using 
the Boltzmann equation 



N 



dt 



V ■ V/w = - ft/iv - Mn. (1) 



where the collision terms P^Ji represent the interaction 
(due to CE and ionization) between the species k and I. 
The interaction rate j3k is formally written as 



/3fc(w,z)= J d WVrel<j{Vrel)fki'W,z) , (2) 

where Vrei — \v ~ w\ and a is the cross section for the 
interaction process. More precisely, (3^ is the rate of 
CE of an ion that becomes a neutral, Pi is the rate of 
CE plus ionization of a neutral due to collisions with 
protons, while (3e is the ionization rate of neutrals due 
to collisions with electrons. Eq. ([Ij is used to calculate 
/at starting from the distribution of charged species (pro- 
tons and electrons) , under the assumption of stationarity 
(dfN/dt = 0).^ 

The dynamics of protons and electrons coupled with 
neutrals can be described very generally through conser- 
vation equations of mass, momentum and energy: 



oz 



Pe)Vi 



= 0, 



(3) 



d-z^^P' 



,)vf + + Pe + F,no,n] = , (4) 



d_ 

dz 



2 7ff 



PeH 



0, (5) 



where F^^ass = mn J d^vv\\fN, Fmom = "iff / d^vvjfN 
and Fen = mHf^ J d^vvn (w| + wl)/Af are respectively the 
fluxes of mass, momentum and energy of neutrals along 
the z direction. Usually the dynamical role of electrons 
is neglected due to their small mass. However, collective 
plasma processes could contribute to equilibrate electron 
and proton temperatures to some level. If this equilibra- 
tion is very efficient, then the electron pressure can no 
longer be neglected and the total gas pressure needs to 
include both proton and electron contributions, namely 
Pg = P, + Pe = -Pz(l + /?), where Piz) = Te/T, is the 
electron to proton temperature ratio and is taken as a 
free parameter. 

The solution of Eqs.(IT])-(l5]) is described in detail in Pa- 
per I. The only difference here is the presence of elec- 
trons, which has been previously neglected because it 
does not affect the shock dynamics, unless ^ Ti. The 
importance of introducing the electron contribution in 
this context comes from the fact that Balmer emission is 
very sensitive to both Tp. and Ti (jHeng fc McCravll2007l : 
Ivan Adelsberg eralll2008[) . Therefore, even very partial 
equilibration can produce observational consequences in 
the line emission, as we show in !21 

The level of electron-ion equilibration is one of the 
open questions in coUisionless shock physics and Ha ob- 
servations present us with a unique tool to investigate 
this aspect. Theoretically the equilibration process is far 
from being understood and anything between total ab- 
sence of interaction and full equilibration has been pro- 
posed (see iRakowskv,,2005L for a review). Moreover the 



electron-ion equilibration level is expected to change be- 
tween upstream and downstream of the shock, because 
the plasma conditions in these two regions are totally 
different in terms of temperature and turbulence proper- 
ties. As a consequence, we distinguish between upstream 
and downstream using two separate parameters, /3„p and 
Pdown, respectively. 

We stress that a crucial assumption of our calculation 
is that new protons produced at position z by CE and 
ionization instantly thermalize with other protons. As 
discussed in Paper I, such an assumption, which is es- 
pecially important in the upstream region, is rather del- 
icate: protons might isotropize and yet not thermalize 
with the rest of protons in the plasma. A dedicated ef- 
fort using numerical particle simulations should be used 
to address this important issue. 

In principle, the dynamics of electrons in the upstream 
region could also be affected by electrons ejected due to 
the ionization of neutrals that have returned upstream. 
This effect can however be neglected for the calculation 
of Balmer lines: in Paper I we showed that the frac- 
tion of returning neutrals is ^ 5%. Let us assume that 
all these neutrals are ionized in the upstream region. 
To a good approximation, stripped electrons can be ex- 
pected to have the same velocity of their parent atoms, 
which, in turn, have a mean velocity of the order of the 
downstream proton thermal speed. Hence the temper- 
ature of stripped electrons is Te.m c± S/SV^jme/fc^ « 
6.2 X 10^(X4h/5000km/s)2K. Now, assuming that the 
newly generated electrons reach equilibrium with pre- 
existing electrons, the final temperature is of the order 
of ("o^o + "■m7i„)/(no -l- Uin) ~ 2.5 X lO^K, where we 
use Vsh = 5000 km/s and an initial ionization fraction 
of 50%. This temperature is too low to be relevant for 
Balmer emission. In fact, electrons can contribute only 
if their velocity is > 2000 km/s, namely if Tg > lO^K. 

3. CALCULATION OF BALMER LINE EMISSION 

Once the distribution functions of hydrogen, protons 
and electrons are known at each position, the calcula- 
tion of line emission is quite straightforward, provided 
that the relevant cross sections are known. Here we con- 
centrate on the Balmer Ha line, which results from the 
hydrogen deexcitation from level 3s and Sd to 2p as well 
as from 3p to 2s. The latter case is complicated by the 
fact that the 2>p level can also decay into Is, producing 
Ly/3 photons. Depending on the optical depth of the 
medium, these photons can either escape the system or 
be reabsorbed by ground-state hydrogen and eventually 
reemitted as Ha photons. In the literature the optically 
thin and optically thick cases are usually labelled as Case 
A and Case B, respectively, and the total Ha production 
rate is written as: 



Rho — Rh{3s 



R 



H{3d) 



B3p,2sRh{3p) 7 (6) 



where Rh{31) is the production rate of hydrogen excited 
at level 3/ and the factor B^p_2s is the fraction of transi- 
tions from 3p to 2s, which is ~ 0.12 in the optically thin 
case (Case A) while it becomes unity in the optically 
thick case (Case B) (jvan Adelsberg etal|[2008l) . 

The conversion efficiency from Ly/3 to Ha photons 
depends on the shock speed, the electron-ion tempera- 
ture ratio and the pre-shock ionization fraction. It was 
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first computed bv lChevalier et all (|1980D . who found that 
in conditions appropriate for many Balmer dominated 
shocks the emission from cold hydrogen is generally op- 
tically thick while the emission from hot hydrogen is close 
to be optically thin. In this work, following that result, 
we adopt B^p^2s = 0.12 (1) for the emission produced by 
hot (cold) hydrogen. 

In order to calculate the Ha emission for both Case A 
and Case B, we need to compute the different production 
rates of hydrogen excited to the sublevels 3s, 3p and id. 
The excited hydrogen is mainly produced by two differ- 
ent processes: collisions with protons and electrons, and 
CE reactions leaving the hydrogen atom in an excited 
state. For the sake of clarity here we neglect further con- 
tributions due to collisions with helium. The production 
rate of H{31) at a fixed position z reads: 



Rh{31){v,z) = J (fwVrelfH{v,z) X 
ft (w, z) cr*°* i)'^ {Vrel ) + fe{w,z) cr*°* g^)'^ {Vrel ) 

+ f d^WVrelMv,z)fH{w,z)<jlf^^''>{Vrel). (7) 



Notice that the hydrogen distribution function fn in- 
cludes only atoms at the ground level Is. In fact, colli- 
sional deexcitation can be neglected because it occurs on 
a typical time scale TcoU ^ iniacoiiVsh)~^ — lO'^sec (for 
Hi = 1 cm~^ and Vsh = 1000 km/s) which is much longer 
than spontaneous deexcitation, whose time scale (from 
state n to state to) is r = h/AE — 4.8 x 10~^^(n^^ - 
TO~^) sec. Hence atoms decay to the ground state before 
undergoing any further collision. 
In Eq. d?]) we use the total excitation and CE cross 

sections for the atomic sublevel 31, al%*^^^^ and a]°e^^^\ 
which take into account the direct excitation processes. 
Is — ^ 3Z, plus contributions coming from atoms excited 
to higher levels with n > 3 that subsequently decay to 
the state 3/. Formally this total cross section can be 
written in the following form: 



tot(3i) 



Br, 



'l',3l '^1.1 



(8) 



where ais^n'i' is the cross section for excitation from 
the ground state to the level n'l' and the Bn'i',ni are 
the cascade matrix elements representing the probability 
that a hydrogen atom excited to state n'l' will make a 
transition to state nl (with n < n') via all cascade routes. 

We treat CE, excitation and ionization between 
electrons, protons and hydrog e n ato ms using the 
cross sectio ns o f iBarnett el al.l (Il990f). lBelki6 et al.l 
(fT99l).|jMK ;v & SmithI (|1993| ). lBalanca et al.l (|1998[ ) a^ 
iHarel et al.l (1199811 ^. For some of these cross se ctions we 
adopt the fi tting functions provided by iHeng fc Sunvae^ 
(|2008| ) and iTsehakhovich et all (|2012l ). At the time of 
writing, the CE cross sections for sublevels with differ- 
ent angular momentum states are known only up to the 
level n = 4. We estimate that neglecting the contribu- 
tion of higher levels entails an error around 5% in the 

^ Many of these cross sections can be found in 
the International Atomic Energy Agency website: 
http : //www-amdis . iaea.org/ALADDIN/ 



calculation of the total <Tce, therefore in the following we 
ignore levels with n > 4. Hence, for the CE process, the 
cross sections reduce to the following: 



tot(3s) 
ce 

tot{3d) 

ce 

*ot(3p) 



(10) 



0'ce,ls~^3d + B4d.3d O'ceAs- 



►4d : 



■ Crce,ls^3p + Bis,3p 0'ce,ls-i.4s 
+ Bid,3p Cce,ls->4d 7 



(11) 

where we use t he v alues of Bi^i 31 as listed by 
IHeng fc SunvaevI (|2008D (see their Table 3). An iden- 
tical approach is adopted to calculate the impact ex- 
citation by protons, restricted to levels 3 and 4. We 
adopt the impa,ct ex citati on cross sections cal c ulated by 
iBalanca et ail (|1998( ) and iTseliakhovich et all (|2012[ ) for 
the sublevels 31 and 4/ respectively. Unfortunately these 
works provide cross sections only for a limited range of 
impact energies, 1-100 keV and 5-80 keV respectively, 
which means that the relative speed between protons 
and hydrogen atoms can be respectively in the range 
[vi,V2] = [438,4377] km/s and [978,3914] km/s. Out- 
side these velocity ranges we estimate the sublevel cross 
sections, g "^^-^, using the total cross sections, (J^^^.-^ as 

taken from lJanev &: SmithI (|1993f ). in the following way: 



(t"L,4«2) 



> = <^u~>iv) ::^\~. (1 + 6„) , (12) 

where the coefficients e„ are chosen in such a way as to 
have J2i '^ni{vi(^2)) = <^n{vi{2)), and their values are of 
the order of few percent. A similar approximation has 
been used also for v < vi. 

Also the impact excitation by electrons is limited to 
sublevels 3/ and 41 and expressions similar to ((9|)- (|lip 

hold for (T*°*/"^/-' as well. In this case we use the cross sec- 

ex{e) 

tions provided by iBrav fc Stelbovicsl (|1995| ). computed 
using the convergent close-coupling method. As for the 
previous case, the error in the total cross section pro- 
duced by excluding higher excited levels is around 5%. 

In order to compute the spatial emissivity profile of the 
Ha emission we need to integrate Eq. (j6]) in the veloc- 
ity directions orthogonal to the line of sight. From the 
observational point of view, most cases refer to shocks 
viewed edgewise because of limb brightening. In such 
cases, assuming a pure plane shock (i.e. neglecting cur- 
vature effects), the edge- wise emissivity profile results 
from the following integration: 



FHaiz,V^) = JJ dv^dvy [Rh(3s){z,v) 



+RH{3d) {z,v) + B3p,2sRh(3p) (z, v)] 

dvj_ 



2 / dv, 

' —00 



[Rh(3s) {z, V^,Vi_) + RH{3d}{z, Vz, V±) 

+ B3p^2sRh{3p){z,Vz,V±)] ,(13) 

where z is the direction of shock propagation, x is 
the direction along the line of sight and y is orthogo- 
nal to the (x, z) plane. In the second equality we use 
vj_ = {v^ + VyY^'^. Starting from Eq. (fT3|) we can obtain 
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more useful integrated quantities, which can be directly 
compared with observations, namely the spatial emis- 
sivity profile, (,Ha{z), the volume-integrated line profile, 
4>Ha{vx), and the total line strength, I Ha, respectively 
defined as 



£,Ha{z) = i dVxFHa{z,Vx) , 
J — oo 

/oo 
dzFHa{z,v^) , and 
-OO 

iHa = dVxdzFHa{z,Va;) ■ (14) 



In the next section we will investigate how these ob- 
servables can be used to test the presence of a neutral- 
induced precursor and to infer the ambient parameters 
(as the ionization fraction, the shock speed and the 
electron-proton equilibration level) . 

4. RESULTS 

In this section we illustrate the main results of our ki- 
netic calculation concerning the Ha emission. It is well 
known that Balmer emission is highly sensitive to plasma 
density, shock velocity, initial degree of ionization and 
electron-ion equilibration level. Our aim is to illustrate 
how to disentangle different effects produced by these 
quantities. The main observable quantities are the spa- 
tial profile, the line profile and the relative intensity of 
the broad and narrow lines. Below we discuss the in- 
terpretation of these quantities and their relation with 
physical parameters, comparing our results with previ- 
ous work. 

For all cases discussed, we fix the temperature at up- 
stream infinity to 10"'K, because larger values are in- 
compatible with the presence of neutral hydrogen, while 
for lower values the results do not change significantly. 
Also the total upstream numerical density is fixed to 
ntot = O.lcm"^, while the plasma ionization fraction is 
assumed to be 50%, unless otherwise specified. We notice 
that the total density can be factorized out in Eqs.(IT])- 
([5]), hence a change in the total density only reflects in 
a change of the length scales of the problem, while all 
other quantities remain unchanged. 

As we already pointed out, the electron-to-proton equi- 
libration level, /3, is expected to change between up- 
stream and downstream of the shock, hence we use two 
independent parameters, (3up and Pdown, respectively. 
We will focus mainly on the two extreme cases of full 
equilibration (FE), where protons and electrons share 
the same temperature everywhere (i.e. Pup = Pdown = 1) 
and no equilibration at all (NE), which corresponds to 
the situation where electrons and protons do not in- 
teract at all. In the latter case the electron tempera- 
ture is equal to lO^K in the entire upstream region and 
Pdown = me/mp in the downstream. Intermediate equi- 
libration cases will be also discussed. 

4.1. Spatial emission 

Fig- [D shows the spatial emissivity profile, ^Ha{z), for 
different shock velocities and for the two extreme as- 
sumptions for the electron-ion equilibration efficiency: 
the upper panel shows the NE case, while the lower panel 
shows the FE case. Interestingly, while in the former 
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Fig. 1. — Spatial emissivity profile of the Ha line for ritot = 
O.lcm^'^, 50% ionization fraction and for different values of shock 
velocity, as specified in each panel. Tv^o extreme cases of electron- 
ion temperature equilibration are shown: NE {upper panel) and 
FE {lower panel), i.e. fiup = l^domn = 1- 



£. 0.3 



-a 0.2 




1000 1500 2000 2500 3000 3500 4000 4500 5000 
Vsh [km/s] 



Fig. 2. — Fraction of the total emissivity coming from the up- 
stream as a function of shock velocity and for three different as- 
sumptions on the degree of electron-proton equilibration upstream 
{Pdo-wn is fixed to 1 for all cases). 
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Fig. 3. — Efficiency of tlie Ha emission (i.e. number of pfiotons 
emitted for eacfi liydrogen atom crossing tlic sfiock) as a func- 
tion of tile sliock velocity and for different degrees of electron-ion 
equilibration efficiency. The tliick lines refer to the cases of non- 
equilibrium (solid), full equilibrium (dotted) and partial equilib- 
rium with Pdown = 0-1 (dot-dashed). For these cases the plasma 
is assumed to be optically thin for Ly/3 photons produced by hot 
hydrogen atoms, but optically thick for photons produced by cold 
hydrogen atoms. The thin lines, instead, show NE cases assuming 
a medium that is either optically thin (Case A) or optically thick 
(Case B) to Ly/J photons. 

case the emission is produced only in the downstream, in 
the latter case a substantial fraction of the Ha emission 
comes fr'om the upstream. This fact is better illustrated 
in Fig. [21 where the fraction of total Ha emission pro- 
duced upstream is plotted as a function of shock velocity 
and for different values of electron-ion equilibration effi- 
ciency upstream (while Pdown is fixed to 1): the upstream 
emission has a peak when Vsh is close to 2500 km/s and 
can reach ^ 40% of the total emission when f3up — 1. 
Remarkably, even reducing /3„p down to 0.05, the up- 
stream emission is still a non-negligible fraction of the 
total, being around 10%. For velocities > 2500 km/s the 
upstream emission decreases monotonically because the 
heating produced by the neutral return flux becomes less 
efficient, as shown in Paper I (see its Fig. 6). 

The Ha emission produced upstream is clearly due to 
collisional excitation of hydrogen by electrons. In fact, 
this process has a threshold for Vrei — 2000 km/s and 
peaks at Vrei — 3000 km/s, hence when the electron tem- 
perature is larger than ~ 1.5 x lO^K the atomic level 
n = 3 is easily excited and Ha emission is produced. On 
the other hand, excitation produced upstream by protons 
is suppressed because, for a given common temperature, 
protons have a thermal speed ^m^jm^ times lower than 
electrons; hence Ha emission induced by proton collisions 
and CE becomes relevant only for Tp > few x lO^K. The 
heating induced by the neutral return flux may lead to 
such ion temperatures upstream, but only very close to 
the shock, on spatial scales that are too small to affect 
the Ha emission. 

From Fig. [1] we can see that Ha emission seems to 
increase with increasing shock velocity. This is a conse- 
quence of the fact that the production efficiency of Ha- 
photons is almost constant at high Vsh- This is shown 
in Fig. |3] where we plot the efficiency of Ha emissivity, 
defined as the number of Ha photons emitted per hydro- 



gen atom crossing the shock, namely e = lual (VshnH.o), 
where uha is the hydrogen number density at upstream 
infinity. In Fig. |3] several cases are shown: the solid and 
dotted thick lines refer to NE and EE cases respectively, 
while the dot-dashed thick line refers to an intermediate 
case (NE upstream and partial equilibrium downstream 
with (idown = 0.1). For aU these three cases we assumed 
that the plasma is optically thin to Ly/3 emission from 
hot hydrogen atoms (Case A) and optically thick to Ly/? 
emission from cold hydrogen (Case B) as discussed at the 
beginning of 521 In order to compare our results with pre- 
vious work, we also plot the cases where the plasma is 
optically thin or optically thick for both the hot and the 
cold hydrogen emission, but only for the NE case. 

These latter two cases can b e compared with the 
estima ted efficiency provided by IChevalier fc Raymond! 
([19781) who give eA = 0.048 and = 0.27, indepen- 
dently of the shock velocity. We note that while es is 
quite close to our finding, especially for high shock speed, 
the value of tA is about a factor two smaller than our re- 
sult. The discrepancy is probably due to the fact that 
the excitation cross sections used at that time were not 
very accurate. 

4.2. Ha lint 'profile 

Let us now analyze the shape of Ha lines. Fig.[31shows 
the volume- integrated line profiles, 4'Ha(vx), for several 
values of the shock velocity, for both NE (upper panel) 
and FE (lower panel) cases. Several comments are in 
order. The first point to note is that (pHa cannot be 
described, in general, by only two Gaussian-like compo- 
nents, as usually assumed in the literature. Besides the 
usual narrow and broad components, we clearly see the 
presence of a third intermediate component whose typi- 
cal width is about few hundreds km/s. This intermedi- 
ate component is a direct consequence of the existence of 
the neutral-induced precursor. In fact, as we showed in 
Paper I. the neutral return flux can heat upstream pro- 
tons up to a temperature ~ 10^ — 10''K. Hence hydrogen 
atoms that undergo CE upstream with these warm pro- 
tons can produce Ha emission with a typical width of 
~ 100 km/s. This picture also suggests that the inter- 
mediate component should have a non-Gaussian profile, 
because it contains contributions from hydrogen popula- 
tions at different locations in the precursor, which have 
different temperatures. 

This physical interpretation of the intermediate com- 
ponent is supported by the fact that its intensity, with 
respect to the narrow line, increases or decreases accord- 
ing to the temperature and length of the neutral pre- 
cursor. For example, in Fig. [5] we show the effect of 
increasing the initial neutral fraction from 1% to 50%. 
We sec that the intermediate component becomes more 
prominent as the neutral fraction increases: this is a con- 
sequence of the fact that the heating induced upstream 
by the neutral return flux increases. A similar behaviour 
occurs when changing the shock speed. In Paper I we 
showed that the upstream temperature has a maximum 
for Vsh — 2000 km/s and decreases for smaller and larger 
speeds. The same happens for the emission of the inter- 
mediate component with respect to the narrow one (see 
Figs, [g and 13). 

The variation of electron-ion equilibration efficiency 
also considerably affects the line profile. In Fig. [6l we 
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Fig. 4. — Volume-integrated Ha emission for different values of 
the shock velocity. The upper and lower panels show the cases NE 
and the FE cases, respectively. The values of total density and 
ionization fraction are the same as in Fig.[T] The small boxes show 
a zoom for small values of Vx , aimed at emphasizing the shape of 
the narrow and intermediate lines. 



0.01 



1e-05 



neutral fraction= 1% 
10% 
50% 



n,ot= 0.1 cm" 
V-h = 2000 km/s 




-3000 -2000 -1000 1000 2000 3000 

Vj, [km/s] 

Fig. 5. — Volume-integrated Ha emission for three different val- 
ues of upstream ionization fraction. 
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Fig. 6. — Volume-integrated Ho emission for different values of 
the downstream electron-ion equilibration efficiency Pdown- The 
upstream electron heating is assumed to be null. 



plot (pHa for a fixed shock speed and for different val- 
ues of Pdown- When j3down increases, a decreasing of the 
broad emission is observed: this occurs because electrons 
contribute to ionize hot hydrogen atoms. Also the width 
of the broad component decreases because a fraction of 
the protons' thermal energy is transferred to electrons, 
hence the proton temperature decreases. The narrow 
component, on the other hand, is only slightly affected 
by variations of Pdown- Its intensity increases when Pdown 
goes from rrie/mp up to 0.01, while for larger values it 
remains constant. In Fig. [5] we only consider the effect 
of electron-ion equilibration downstream, while the elec- 
tron temperature upstream is taken constant and equal 
to lO^'K. The effect of electron heating upstream can be 
appreciated by looking at Fig. [7l where we plot sepa- 
rately the volume-integrated emission from upstream and 
downstream assuming FE downstream, but distinguish- 
ing the NE and the FE cases upstream. The FE case 
shows that the total upstream emission is comparable to 
the downstream one, but has a very different line profile, 
which strongly departs from a Gaussian shape. More- 
over, no broad line comes from the upstream. 

In order to perform a more quantitative study of the 
Ha emission with the aim of extracting independent in- 
formation from the three components, we decide to fit 
the whole line profile using three Gaussian curves. Some 
examples of best fit profiles are shown in Figs. [8] and [9] for 
the NE and the FE cases, respectively. The first point to 
notice is that the shape of both the broad and intermedi- 
ate components slightly departs, in general, from a per- 
fect Gaussian profile. This is especially true for the NE 
case and for a shock speed value below 2500 km/s (see up- 
per panels of Fig. [8]) . On the other hand when the shock 
speed and/or the electron-ion equilibration level increase, 
the quality of the fit improves noticeably. We note that 
the deviation of the broad component f rom a pure Gaus- 
sian profile was already pointed out by iHeng &: McCravl 
( 2007) . 

Using the 3-Gaussian fit, we extract the FWHM of all 
three lines, which provides several pieces of information. 
The first remarkable result is that the width of the nar- 
row component docs not change significantly varying the 
shock speed and the initial ionization fraction. Its value 
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Fig. 7. — Volume-integrated Ha emission from downstream (solid 
lines) and upstream (dashed lines) for two extreme cases of up- 
stream electron heating: thin lines refer to the case with constant 
upstream electron temperature, Te = 10'*K, while thick lines are 
computed for complete equilibration jSup = 1 . The value of down- 
stream equilibration is fixed to Pdown = 1 for both cases. 
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Fig. 8. — Examples of fits with three Gaussian curves for different 
values of the shock velocity, in the NE case. The dotted lines 
represent the emission computed with our kinetic model, while the 
solid line is the best fit. The quality of the fit increases for larger 
values olVsh- Notice that all profiles are normalized such that the 
peak value is equal to 1. 



is always ~ 21 km/s, which corresponds to a population 
of atoms with a temperature of 10^ K. This result im- 
plies that the neutral-induced precursor does not affect 
at all the narrow line width, which is only determined 
by the hydrogen temperature at upstream infinity. This 
is a consequence of the fact that the precursor length is 
smaller than the CE interaction length of cold hydrogen 
atoms in the upstream, irrespective of Vsh and other am- 
bient parameters. This result is particularly important 
because it demonstrates that the neutral precursor can- 
not be responsible for the anomalous narrow-line com- 
ponent detected from several SNR shocks, which have a 
FWHM as large as -30-50 km/s. 

Concerning the broad and the intermediate compo- 
nents, their FWHM are shown in the upper and lower 
panel of Fig. [TUl respectively. As usual we plot the re- 
sults for the NE and FE cases plus some intermediate 
cases for the electron-ion equilibration level. Our results 




Fig. 9. — The same as Fig.|8]but for the FE case. 



for the FWHM of the broad compone nt are in good agree - 
ment with previous calculations (^e.g. lSmith et al.iri991|) . 
The only exception concerns the NE case, which departs 
from the general trend for Vsh < 2500 km/s. This be- 
havior does not have a direct physical meaning and, as 
already noticed, is rather due to the particularly bad 
quality of the fit in this region of the parameter space. 
We remark that the quality of the fit rapidly improves 
for larger value of /3 and the broad line width resulting 
from the fit becomes very close to the actual width when 

f^down ~ 0.01. 



As first pointed out by I Chevalier et alj (|198CI| ). the 
FWHM of the broad line is a direct measurement of the 
proton temperature downstream. As a consequence it 
only depends on the values of Vsh and Pdown ■ This result 
can be easily understood using a plane parallel shock 
model for a totally ionized plasma, which gives a pro- 
ton temperature equal to Ti = 'impV'^j^/\(S{l + Pdowajkn, 
where nip is the proton mass and Ub the Boltzmann con- 
stant. As we showed in Paper I, this result still remains 
a good approximation when the plasma is only partially 
ionized and the neutral return flux is taken into account. 
On the other hand, a deviation of the proton temperature 
from this estimate can be induced by the presence of he- 
lium. In fact, immediately downstream of the shock, he- 
lium nuclei thermalize at a temperature mfje/'^p times 
larger than the protons' one. If helium and protons ther- 
malize at the same temperature on a length scale smaller 
than the excitation length scale, the mean temperature of 
hot hydrogen produced by CE with hot protons is larger 
than the prediction without helium. As a consequence 
the FWHM is expected to be larger. Indeed, the presence 
of heli um was taken into account bv lvan Adelsberg et al.l 
(j200^ , which found for the broad component a FWHM 
about 15-18% larger than the one found in our calcula- 
tions. 

Let us now consider the width of the intermediate com- 
ponent (lower panel of Fig. [TO]). In this case, for 1000 < 
Vsh < 5000 km/s, the FWHM ranges between 100 and 
300 km/s. A peak is present for 2000 < Vsh < 3500 
km/s, depending on the level of electron-ion equilibra- 
tion. Once again, we notice that for Vsh < 1500 km/s 
the FWHM obtained from the fit procedure is not well 
determined in that the intermediate component departs 
from a pure Gaussian shape. Moreover, the emission due 
to the broad component becomes much larger than the 
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contribution of the intermediate one, which, in turn, be- 
comes less distinguishable. 

Observational evidences, compatible to what we call 
here intermediate component^ have been reported in sev- 
eral works, even if such evidences have never been related 
to the neutral-induced precursor. The most interesting 
case is the Ha line pr ofile detected frorn the " Knot g" 
of the Tycho's SNR bv lGhavamian et al.l (|2000f ). There, 
an observation of Ha emission performed with high spec- 
tral resolution suggests the presence of two superimposed 
lines: a narrow one, with a FWHM of 44 km/s, plus 
a second, less pronounced line, whose FWHM is ^ 150 
km/s. Such a width is fully compatible with the FWHM 
of the intermediate component resulting from our cal- 
culations. On the other hand it is important to stress 
that the Knot 3 is a complex region where density vari- 
ations of the ISM produce a non-spherical shock, hence 
the observed line profile could also result from projection 
effects. For this reason observations with better spatial 
and spectral resolution are required in order to disentan- 
gle geometrical effects from physical effects. The Balmer 
emission detected from Tycho is not an isolated case. 
Several spectra observed from Balmer-dominated shocks 
show evidence of narrow Ha lines with non-Gaussian 
"wings" (see e.g. iSmith et al.l 11994 ). We suggest that 
such wings could be the signature of the intermediate 
component. 

In spite of this interesting connection, in the present 
work we avoid performing a detailed comparison between 
theoretical predictions and observations. The reason is 
that although the calculations presented here represent a 
considerable step forward in the description of collision- 
less shocks in partially ionized media, they still remain 
incomplete. Aside from some minor complications that 
need to be taken into account (like the presence of he- 
lium or projection effects arising from deviations from 
plane geometry), a major role in determining the shape 
of the Ha line is played by the presence of CRs, which 
are not yet included here. In fact, it is widely accepted 
that efficient CR acceleration changes the global shock 
structure, generating an upstream precursor that is, to 
some extent, similar to the neutral- induced precursor but 
develops on a different length scale. Effects induced by 
CRs could be especially relevant for the Tycho's SNR, 
which has been suggested to accelerate CRs efficiently 
(jMorlino fc Capriolill20T2l ). In passing we notice that the 
efficient acceleration of CRs is the most plausible expla- 
nation of the relatively wide narrow Balmer line found 
in Tycho (FWHM of ~ 44 km/s). 

It is worth noticing that alternative explanations for 
the non-Gaussian wings have also been proposed. For 
example, Jlaymond et al. (2008) suggest that deviations 
from the Gaussian profile could be the result of a non- 
Maxwellian proton distribution downstream. In fact, 
neutral atoms that become ionized could settle into a 
bi-spherical distribution (similar to that of pickup ions 
in the solar wind) that would then introduce a non- 
Gaussian contribution to the line core. We recall that 
in the present work, as well as in Paper I, we do not 
include such an effect, but we assume that, soon after 
being ionized, atoms thermalize with the rest of ions. 

4.3. Ha emission from the upstream 
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Fig. 10. — FWHM for broad (upper panel) and intermediate com- 
ponents {lower panel) resulting from the the fit of total Ha emission 
with 3-Gaussian distributions. 



This section is devoted to highlight some features of the 
Ha emission from the upstream region. This is a crucial 
observable in order to understand the effects produced by 
CR acceleration, as we will show in a forthcoming paper. 
In the near future, observations could reach good enough 
spatial and spectral resolution so as to provide detailed 
spectra at different distances from the shock position, 
which makes the exercise of analyzing the details of the 
emission from the upstream especially interesting. 

As we already showed in Fig. [2] the upstream region 
could radiate up to 40% of the total Balmer emission 
in the case of full electron-ion temperature equilibration. 
On the other hand if /3„p ^ 0.01 the upstream emission 
drops below 1% of the total. The line profile in the up- 
stream emission is quite different from the one produced 
in the downstream region. In fact only the narrow and 
the intermediate lines are present. This is clearly shown 
in Fig. [H] where the upstream line profile at different 
distances from the shock is shown for the cases of FE 
and partial equilibration and for Vsh = 2000 and 4000 
km/s. We chose the following distances: z = lx,2x 
and 5 x L^fp, where L^fp = l/(a'ce«tot) ~ 10^*' cm, 
and the charge-exchange cross section is approximated as 
fJce ^ 10~^^ cm^. As we move far away from the shock, 
the FWHM of the intermediate line decreases as a conse- 
quence of a reduction of the temperature in the precursor, 
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Fig. 12. — Spatial cmissivity profile of the Ha line in the up- 
stream. The upper panel shows the case of partial electron-ion 
equilibration with /3up = fidown = 0.1 while the lower panel shows 
the FE case. Dashed and dot-dashed lines show the contribution 
of narrow and intermediate lines, respectively, while the solid line 
is the total emission. The line thickness (and color) distinguishes 
the different shock velocities: 1000 (thin-red), 2000 (middle-green) 
and 4000 km/s (thick-blue). 



Fig. 11. — Emissivity profile at three different locations in the up- 
stream for partial equilibration, fiup = Pdown = 0.1 {upper panel), 
and FE, /3up = Pdown = 1 {lower panel). Thick (red) lines show the 
ease for V^h = 2000 km/s while thin (blue) lines are for Vah = 4000 
km/s. 

while the narrow hne has always the same FWHM. The 
total emission falls down at a distance of ^few Lmfp- This 
distance corresponds to the position where the electron 
temperature falls below the threshold of « 1.5 x lO^K 
discussed in the second paragraph of § 14.11 This point 
moves further from the shock for larger values of (3up as 
can be clearly seen in Fig. [T^l where we plot the inte- 
grated line emission as a function of the position only 
in the upstream, distinguishing the contribution of the 
narrow and intermediate line and for different values of 
the shock velocity. When we have FE the contribution of 
the intermediate line is always smaller than that of the 
narrow line, but for lower equilibration levels the ratio 
of intermediate over narrow emission increases and for 
Pup = 0.1 the two lines contribute at roughly the same 
level. These findings are summarized in Fig. 1131 

4.4. Simulated observations 

Although the three-Gaussian fit presented in Section 
§ 14.21 catches the essence of the expected distorsions in 
the Balmer lines, it does not reflect the whole complex- 
ity involved in fitting observed Ha line profiles. Actual 
data are affected by a number of instrumental and sta- 
tistical issues, the most obvious of which (and the only 
ones that we will investigate here) are: the limited spec- 
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Fig. 13. — Ratio of the intermediate to narrow Ha emission at 
different locations in the upstream, as a function of the shock ve- 
locity. Different line shapes correspond to different locations, as 
indicated in the label. Thick (blue) lines show the partial equili- 
bration case, while thin (red) lines show the FE ease. The region 
of the parameter space where the Ha emission drops at a level of 
10""* times the peak value are not shown in the plot. For exam- 
ple, for the case of partial equilibration at a distance of SLj^fp the 
emission is always below this threshold. 
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Fig. 14. — Results of a 2-Gaussians fit to a simulated line profile, 
obtained by downgrading our Model with a typical instrumental 
resolution and photon noise, and adding Poisson noise (see text for 
details). The upper panel shows the simulated data (soUd line), 
together with the best-fit narrow line (dashed line) and broad line 
(dotted line) components; while the lower panel shows the residuals 
after subtraction of the best-fit solution. 



tral resolution of the instrument, the Poisson noise of 
the line photons themselves, and any additional photon 
noise, either due to the astronomical or to the instru- 
mental background. 

For the simulations presented here we have used the 
line profile computed for Vsh = 3000 km/s (and with an 
ion fraction of 50%), Pdown = "m^jm^^ re(up)= 10^ K. 
This situation is the one plotted as a dotted line in the 
upper panel of Fig. |4l 

As for the observational parameters, we will adopt a 
published observation o f an Ha line profile in SN 1006 
(jGhavamian et all 120021) as the reference for the instru- 
mental parameters as well as for the flux levels. In that 
observation, the instrumental resolution is 4.5 A, corre- 
sponding to 205 km/s, while the dispersion per pixel is 
0.27 times the resolution; the total number of photons 
measured in the line is about 10^ (for a 140 min integra- 
tion time, with a 4-m class telescope), while the back- 
ground noise level is about 3000 photons/ A: these photon 
numbers are for the whol e spectrograph slit and in tegra- 
tion time, as specified by iGhavamian et al.l ()2002|) . 

Fig. [T3] shows the results of a 2-Gaussian fit to simu- 
lated data, obtained combining the model with the in- 
strumental parameters mentioned above. As shown by 
the residuals, in this case the quality of the observation 
is not sufficient to investigate details, beyond the mere 
separation of a narrow and a broad component. 

We have then explored several combinations of the ob- 
servational parameters, focusing on the spectral resolu- 
tion (expressed in terms of velocity resolution, Vres) and 



on the photon statistics (expressed in terms of the to- 
tal number of photons in the line, Nphot)- As for the 
instrumental dispersion per pixel, we have kept the ra- 
tio of 0.27 times the resolution, as in IGhavamian et all 
(l200l . while we have usually adopted the background 
noise level given above. We have also tried with a much 
lower background noise level, but, of course, even in this 
case, the noise component associated with the photons 
of the line itself cannot be eliminated. 

We have used a grid of simulations to investigate, on 
the log(iVp/iot)-Wres parameter plane, the behavior of 2- 
Gaussian and 3-Gaussian fits. In both cases we have 
chosen a grid of 25 x 25 points, suitably positioned in the 
parameter plane. In order to minimize the "noisy" look 
in the figure (a natural consequence of the fact that each 
simulation includes random noise), we have performed a 
large number (100) of simulations for each point; out of 
them, we have discarded the 5 cases with the highest 
value as well as the 5 ones with the lowest value, and 
we have then taken the average x^ of the remaining ones. 
The results are shown in Fig. [151 

As a result, for our Model, we may see that, for a line 
photon statistics of about 10^ photons, a spectral resolu- 
tion better than about 180 km/s is required to show (at 
a 3-(T confidence level) that a 2-Gaussian fit is not ade- 
quate; while a resolution better than about 70 km/s is 
required to challenge the 3-Gaussian fit. Even if having 
more photons does matter, in general the photon statis- 
tics does not seem to be a parameter as crucial as the 
spectral resolution. Of course, in order to resolve the 
narrowest spectral component a considerably higher res- 
olution is required; otherwise, it will be detected only as 
an "unresolved spectral component" . 

4.5. Line intensity ratios 

Another observable that may be useful in order to con- 
strain shock parameters is the ratio between the intensi- 
ties of broad and narrow components, h/ In- This infor- 
mation is usually used in combination with the FWHM 
of the broad component, in order to infer simultaneously 
both Vs ft and the le vel of electron- ion equilibration down- 
stream (lHendl2009n . The presence of the neutral-induced 
precursor complicates a bit this exercise because the up- 
stream equilibration also plays a role, as we show below. 

At this point we need to comment on an observational 
caveat. When Balmer emission is observed with a high 
spectral resolution, in order to resolve the narrow com- 
ponent, usually the broad component i s not detected due 
to th e high spectral dispersion (se e.g. IGhavamian et all 
l2000f) . In order to measure the intensity of the broad 
component, which allows one to estimate the It/In inten- 
sity ratio, observations must be performed instead with 
a lower spectral resolution, typically equivalent to a ve- 
locity resolution Av ~ 100 — 300 km/s; but this does not 
allow to resolve simultaneously all three components, be- 
cause at such resolution the intermediate component can- 
not be distinguished from the narrow one, as they appear 
as a single line. As a consequence, in order to provide 
an estimate of lb/ In, we first convolve the line profile 
obtained by our kinetic calculation, with the typical in- 
strumental resolution, Av. Then we fit the convolved 
emission with a two-Gaussian profile, evaluating both 
and /„ from the fitting curves. We choose Av = 150 
km/s, which corresponds to a wavelength resolution of 



12 



Morlino et al. 




8.0 8.2 8.4 8.6 8.8 9.0 9.2 9.4 
Log(N photons) 



160 
140 
„ 120 
I 100 



X 

^ 80 



60 



40 











' .' 






c ^ , 






( ) ; 












/ ' ' 






J t ; 












r .' 
r ' ' 












} / 



























































8.6 8.8 9.0 9.2 9.4 9.6 9.8 10.0 

Log(N photons) 



Fig. 15. — Behaviour of the 2-Gaussians (upper panel) and 3- 
Gaussians (lower panel) fits on simulated data obtained start- 
ing from our Model. We have investigated the parameter plane 
logj^Q (A'j,;ioi)-Vres- The three contours in each panel, drawn re- 
spectively as a solid, a dashed and a dotted line, indicate the 1-cr, 
2-(T and 3-o" confidence level. The contours are obtained by inter- 
polating over a 40 X 40 data grid, each point of which is an average 
value of 50 independent simulations. These plots show that only 
on the lower right side of the panels the quality of the observations 
is good enough to allow outlining deviations from, respectively, a 
2- or 3-Gaussians profile. 



AA = 3.3 A. 

The results are shown in Fig. [1^1 where different lines 
represent different assumptions for the electron-ion equi- 
libration level. The qualitative beha viour is similar to 
what predicted by pre vious studies ([Smith et al.l 119911 : 
IHeng fc McCra^ [20071 ) ■ even if some differences can be 
noted. Our results are, in general, a factor 2-3 larger 
than what predicted by iSmith et al.l (I1991D (see their 
Fig. 8). Their equilibrated case (which corresponds to 
our dot-dashed line) peaks at V^h = 2000 km/s and is 
Ihjln — 2 while our curve peaks at Vgh — 1500 km/s 
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Fig. 16. — Ratio of the broad to narrow Ho emission versus shock 
velocity after the convolution with experimental velocity resolution 
of Av = 150 km/s. Different lines represent different cases of 
electron-ion equilibration efficiency. Thick curves all assume non- 
equilibrium upstream, while the thin curve shows the FE case. 



with lb/ In — 5. In fact, it is rather difficult to per- 
form a close comparison between the ratios computed 
by different authors, because of substantial differences 
in the model assumptions, in the methods used, in the 
assumed chemical composition, and even in the cross sec- 
tions adopted for the various processes: therefore we take 
the above level of differences as acceptable. 

From the observational point of view, in all the SNRs 
for which /(,//„ has been measured, an intensity ratio 
above 1.2 has never bee n seen, while in most cases it 
falls below unity (see e.g. IHeng fc McCravl 120071 ) . If we 
assume no equilibration upstream, this result suggests an 
intermediate value for the clcctron-ion equilibration effi- 
ciency downstream. On the other hand we also see that 
the FE model (both upstream and downstream) predicts 
an intensity ratio < 1 for all shock velocities considered. 
Unfortunately, for given values of lb /In and Vsh there is 
a degeneracy for the values of /3„p and Pdown- 

Moreover, the trend of h/In with respect to the 
electron- ion equilibra ti on is not monotonic. As first no- 
ticed by ISmith et all ([1991D . NE and FE cases do not 
represent the extreme values of the intensity ratio. We 
see, in fact, that for intermediate values of (3 down, h/In 
drops below the equilibrated case Pdown = 1- 

5. CONCLUSIONS 

In this paper we computed the Ha emission produced 
by a coUisionlcss shock that propagates into a partially 
ionized medium. In order to do this, wc first derived the 
evolution of the various species across the shock, by using 
the kinetic model developed in Paper I. This model ap- 
plies to plane-parallel non-radiative shocks in the steady 
state where ions are treated as a fluid, while neutral par- 
ticles are described using the full 3D velocity distribution 
function. On top of this solution we then computed the 
Balmer emission produced by collisional excitation of hy- 
drogen atoms with both ions and electrons, as well as by 
CE events leading to neutrals in excited states. Results 
for the spatial emission and for the line profile of Ha are 
presented for a shock seen edge-wise, varying the shock 
speed, the initial ionization fraction and the electron-ion 
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equilibration level. 

According to the traditional picture 

(|Chevalier fc Ravmon"dl I1978D . the Ha profile de- 
tected from Balmer dominated shocks usually consists 
of two components: a narrow one, whose width reflects 
the temperature of the upstream medium, and a broad 
one, due to neutrals that have undergone CE with hot 
protons in the downstream region. 

This picture is however an oversimplification of what 
happens in reality, mainly for two reasons: i) it assumes 
that neutrals can be described as a fiuid, namely with 
Maxwellian velocity distributions; ii) it does not take 
into account the effect induced by the neutral precur- 
sor. In fact, the latter point is a natural consequence of 
the former one: already in Paper I we showed that, for a 
wide range of shock velocities, a fraction of the hot neu- 
trals produced downstream can recross the shock toward 
upstream, giving rise to a neutral return fiux, and that 
the interaction of these neutrals with the upstream ions 
produces a precursor region where the incoming plasma 
is heated and slowed down. In this paper we have shown 
that the neutral-induced precursor is responsible for the 
production of a new Ha line component, whose width 
is intermediate between the narrow and the broad lines, 
being around few hundreds km/s for a shock speed of 
a few thousands km/s. This intermediate line is due to 
cold hydrogen atoms that have undergone CE with warm 
protons in the neutral-induced precursor, hence its width 
reflects the mean temperature of the precursor. In addi- 
tion, we found that the profllcs of the intermediate and 
broad-line components may deviate from pure Gaussians, 
and that these deviations could be detected by carrying 
out observations of suitably high quality. 

A natural question to ask is whether the heating pro- 
duced in the neutral precursor could explain the anoma- 
lous width of narrow lines, which has already been ob- 
served in so me Balmer filaments associated with s everal 
SNR shocks (|Smith et al.lll994HHester"eraII[T99l . Our 
results show that this is not the case: the bulk of in- 
coming neutrals does not interact with ions in the neu- 
tral precursor, because its extent, which corresponds to 
the interaction length of the returning neutrals, is much 
smaller than the CE length of the incoming neutrals. In- 
stead, as we already discussed, the fraction of incoming 
neutrals interacting with ions in the precursor will give 
rise to the intermediate Ha line. Therefore, we conclude 
that other mechanisms, like for instance a CR-induced 
precursor, should be invoked to explain an anomalous 
width of the narrow- line component. 

Remarkably, some observations point towards the ex- 
istence of intermediate lines compatible with our predic- 
tions: narrow lines detected from Tycho and from other 
SNRs show non-Gaussian wings which could be explained 
with the existence of a third line component. At the mo- 
ment this result must be taken with care because projec- 
tion effects could also be responsible for the observed line 
profiles. Better spatial and spectral resolution are needed 
to disentangle these effects. Unfortunately, the majority 
of these observations do not have the spectral resolution 
required to carry out a satisfactory 3-component fit. In 
order to estimate the experimental requirements neces- 
sary to identify the intermediate line, we compared our 
theoretical Ha profile with simulated observations, which 



take into account both instrumental resolution and Pois- 
son noise of the line photons. As a result, for a typical 
line photon statistics of about 10^ photons, a spectral res- 
olution better than ^ 70 km/s is required to separately 
identify all three components. 

The presence of the intermediate line component may 
also affect the evaluation of the broad to narrow line 
intensity ratio, h/In- This quantity is usually used to- 
gether with the broad line width, in order to estimate 
the level of electron-ion temperature equilibration in the 
post-shock region. Evaluation of lb/ In is usually done 
from observations with resolution > 100 km/s, necessary 
to detect both the narrow and the broad lines. 

This implies that the intermediate line is not resolved 
and that its emission contributes to the observed nar- 
row line intensity. We have included this effect in 
the evaluation of It/ In, and have shown how this ra- 
tio changes varying the electron-ion equilibration down- 
stream. As already pointed out by several authors, if ef- 
ficient electron-ion equilibration occurs downstream, the 
width of the broad line is reduced because a fraction of 
the kinetic energy of incoming protons is transferred to 
electrons. 

In addition to the effects produced by electron-ion equi- 
libration downstream, we investigated what happens if 
equilibration also occurs in the precursor region. Inter- 
estingly, we showed that, increasing the efficiency of equi- 
libration beyond a few percent, electrons can coUisionally 
excite hydrogen atoms, giving rise to Balmer emission 
also from the precursor region. For Vsh ~ 2500 km/s, if 
equilibration is complete (Tg = Tp), the emission from 
the precursor can contribute up to ~ 40% of the total 
Ha emission. This result can be instrumen t al to explain 
the results recently published bv lLee et al.l (|2010l ). They 
observe the Eastern limb of Tycho's SNR, finding a grad- 
ual increase of Ha intensity just ahead of the shock front, 
which they interpret as emission from a thin shock pre- 
cursor. They estimate that the precursor emission may 
contribute up to 30%-40% of the total narrow component 
emission and suggest that the precursor is likely due to 
CRs. In light of our results it is clear that a correct in- 
terpretation of the pre-shock Ha emission requires also 
the evaluation of the emission produced by the neutral 
induced precursor. 

In a forthcoming paper, currently in preparation, we 
will describe the theory of coUisionless shocks in partially 
ionized media in the presence of accelerated particles that 
exert a pressure on the incoming ions. In other words, 
we will generalize the non-linear theory of particle ac- 
celeration in coUisionless shocks to include the neutral 
return flux discussed in Paper I. In the same paper we 
will calculate the shape of the Balmer lines as they are 
affected by accelerated particles, and show how to use 
the widths of the narrow, intermediate and broad com- 
ponents of the Balmer line as tools to measure the CR 
acceleration efficiency in SNRs. 
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